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Abstract 

Master equations describe the quantum dynamics of open systems interacting with an 
environment. They play an increasingly important role in understanding the emergence 
of semiclassical behavior and the generation of entropy, both being related to quantum de- 
coherence. Presently we derive the exact master equation for a homogeneous scalar Higgs 
or inflaton like field coupled to an environment field. It is represented here by an infi- 
nite set of harmonic oscillators with a given spectral density. Our nonperturbative result 
follows directly from the path integral representation of the density matrix propagator. 
Applications and generalizations are discussed. 



1 Introduction 

Often the physical system under consideration cannot be studied in isolation. In such cases 

the interaction with the environment has to be taken into account from the outset, despite 

the fact that one would like to eliminate these 'uninteresting' degrees of freedom from the 
*A preliminary version will be published in Proceedings of Hadron Physics VI, E. Ferreira et al., eds. (World 
Scientific, Singapore, 1999). 



1 



The Functional Derivation of Master Equations 



2 



dynamical description of the proper open subsystem one is really interested in. Prominent 
classical examples are the motion of a massive body under the influence of friction or the 
Brownian motion, where the idealization of an isolated system clearly misses the relevant effects. 
It is not hard to imagine corresponding quantum mechanical situations, most notably the effect 
of Ohmic resistance to the electron current in a wire, and similar other dissipative systems with 
an interplay of transport and relaxation phenomena. The results of research in this field of 
the so-called quantum Brownian motion have been reviewed, e.g., in Ref. PJ, where numerous 
further references can be found. 

Recently open systems have been studied also in the context of quantum field theories, 
with applications ranging from high-energy physics to cosmology. The underlying theme herein 
is the study of quantum decoherence, i.e. the emergence of semiclassical properties in com- 
plex quantum mechanical systems. In terms of the density matrix of the open subsystem one 
would like to understand how a more or less pure initial state evolves into an impure mixture. 
Clearly this question is intimately related to the problem of entropy generation and possibly 
thermalization in quantum systems. The references illustrate the different physical aspects of 
the environment induced decoherence 0, EjL 01 p|, which also may solve the (in)famous problem 
of the 'wave function collapse' || 0. 

In any case, let us recall the general structure of all master equations as expressed in terms of 
density matrices. - The overall density matrix g describing the subsystem plus its environment 
can be expanded in a complete set of states \l/ n , which obey the time dependent Schrodinger 
equation, g = X) n c n |\l/ n )(\l/ n | . Then, employing the Schrodinger equation, the equation of 
motion for the density matrix follows immediately: id t g = [Hse, q] , where H S e denotes the 
total system plus environment Hamiltonian; our units are such that h = c = ks = 1 ■ The 
aim of all derivations of master equations is to reduce the equation for g to an equation for 
the subsystem density matrix p , which is obtained by eliminating the environment degrees of 
freedom, p = TieQ ■ The result, generally, is of the form 



where H is the subsystem Hamiltonian and the remaining term describes the dissipative and 
diffusive effects of the integrated out environment; as we shall see in the following, it may 




(1) 
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have a complicated time dependent (memory) structure possibly involving derivatives of the 
subsystem density matrix p . 

In order to proceed, we consider a simple model motivated by studies of the inflationary 
universe || P . The semiclassical homogeneous field equation for the scalar 'inflaton' $ is 

• • dV 

$ + + r$ = -— , (2) 

(9$ 

where H = R/R is the Hubble 'constant' and the corresponding friction term here describes 
the damping due to the expansion of the universe; similarly, the term involving the decay rate V 
describes the essentially dissipative coupling to other particles; the r.h.s. contains the ordinary 
force term due to the selfinteraction, e.g. V(&) = — |/i 2 $ 2 + |A$ 4 . 

Whereas the fate of the inflaton field driving the early universe expansion has been much 
studied including its quantum as well as thermal fluctuations, we would like to draw attention to 
another aspect of Eq. (|2|). It could arise from a typical semiclassical saddle-point approximation 
to a path integral describing the evolution of $ in a higher-dimensional space-time environment 
[[]. |l(J. In this case, the interesting long-time behavior of the inflaton field will be modified by 
the influence of quantum fluctuations in the higher-dimensional background space. 

In the following, leaving the wider scope for now, we consider as our starting point the 
classical action 

^=/dt|i$ 2 -v($)+f:^[0 n 2 -^(0 n -$) 2 ]} , (3) 

i.e. coupling the homogeneous scalar field to an environment of iV — > oo oscillators, which is 
equivalent to a single radiation field [TJ, |3|. This Feynman- Vernon or Caldeira-Leggett type of 
model represents a large variety of different physical situations depending on the distribution 
of oscillator masses (dimensionless) and frequencies, and the potential, m n , u n , and V(§), 
respectively. 

Similar models have been studied in numerous special cases before, see e.g. Refs. [p], [|, ||. In 
particular, the effective ^-equation obtained in the saddle-point approximation after integrating 
out the environment degrees of freedom exactly (cf. Sec. 2) assumes the form of Eq. (§) for the 
so-called Ohmic environment cf. Sec. 3. 

In distinction to most earlier work, we do not attempt to solve the quantum mechanical 
time evolution problem for the density submatrix of the $-field in a special case. Rather our 



The Functional Derivation of Master Equations 



4 



aim is to derive its master equation employing the path integral representation; see Ref. 
for comparison, where a master equation was derived using the canonical quantization of a 
somewhat different model. We believe that the nonperturbative method presented here lends 
itself to generalizations for more genuine field theoretic problems, i.e. involving propagating 
fields |J, provided the relevant influence functional (cf. Sec. 2) can be calculated. 



2 From the Influence Functional to the Master Equation 

To begin with, since the above action, Eq. (|3|), is quadratic in the oscillator coordinates and 

velocities 4> n , <p n , they can be eliminated exactly from the path integral (s) describing the time 

evolution of the complete system In the following we neglect initial state correlations 

between the proper subsystem and the environment degrees of freedom for simplicity; their 

incorporation does not add a new element to our derivation. 

It turns out that the result of the Gaussian integrations can be specified completely in terms 

of the thermal noise and the dissipation kernels, v and n, respectively, 

f°° (3lu f°° d 

vis) = dtu I(co) coth( — ) cos(o;s) , r](s) = — / dto I(uo) sin(cjs) = — n(s) , (4) 
Jo 2 ' Jo as 

where (3 = T -1 denotes the inverse temperature. The spectral density here is defined by 

1 N 

J M = 9 ^m n ul5{u-u n ) . (5) 

/ n=l 

It characterizes the particular environment. Assuming a quasi-continuous distribution of envi- 
ronmental oscillators, spectral densities I{oS) oc uo k (for sufficiently small oS) with k > have 
often been considered before, cf. |], 0, 3, and have been called sub-Ohmic (k < 1), Ohmic 
(k = 1), and supra-Ohmic (k > 1), respectively. The electron with its radiation field, for 
example, amounts to a particle in a supra-Ohmic environment with k — 3 ||. 

In this way, the following representation for the $-'field' density matrix at time t in terms 
of the initial one (t = 0) is obtained: 

P ($, t) = j d^d^ j($, t- $ i5 % o) P ($i, o) , (6) 

with the propagator determined by two remaining path integrals Jl|, 

J($,$ / ,t;$ l ,$:,0) = ^ /"*'*'* VqVq' exp{i(S[q] - S[q'])}ex V {-F[q,q',t}} , (7) 
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and with the indicated boundary conditions, q(0) = q{t) = and for q' analogously. Here 
Z is a (generally time- dependent) normalization factor, such that the total probability to find 
the subsystem is preserved (see the first of Refs. ||] for a detailed discussion of general properties 
of open quantum systems, the subsystem density matrix in particular, as well as an explicit 
calculation); we omit Z in the following, since it can always be restored by properly normalizing 
p in the very end. Furthermore, the action S is the part of Sse, Eq. ([|), which only depends 
on the field $ , and the Feynman- Vernon influence functional T incorporates the effect of the 
integrated out environment on the subsystem. In the present case it is given by 

F[q,<l',t] = I ds I du[q{s) - q'(s)]{u(s -u)[q(u) - q'{u)] + ir](s -u)[q(u) +q'(u)]} 
Jo Jo 

+ifj(0) r t ds[q 2 (s)-q'\s)] . (8) 
Jo 

Clearly, the spectral distribution has to vanish sufficiently fast in the infrared, in order to 
avoid a divergence of f)(0), cf. Eqs. (f|). We observe that the last two terms in Eq. (§) simply 
renormalize the quadratic parts of the potentials in S[q] — S[q f ] in Eq. (J?]); therefore, we omit 
these terms in the following and add them again in the end, see Eqs. fll4|) . 

Having set the stage with Eqs. — (Rh , we derive the master equation by suitably evaluating 
the time derivative d t p = lim e ^o[p(t+e)— p(t)]/e. We calculate to 0(e) the r.h.s. of p($, $',£) + 
p($, t)e ~ p($, t + e) by going through the following steps: 

i. We split the path integrals going from to t + e at time t, cf. Eqs. by inserting the 
intermediate 'points' z, z' and integrating - corresponding to inserting a complete set of states 
in the formula exp{iH(t + e)} = J dz exp{iHe}\z) (z\ exp{iHt} . 

ii. We parametrize the (infinitesimal) propagation from t to t + e by straight line paths, 

q(s) =z + ^^(s-t) , 'q( s ) = ^—^, t<s<t + e, (9) 
e e 

and analogously q' ; this parametrization becomes exact as e —>■ . 

iii. We introduce x = z ~ ^ > an d x' correspondingly, and replace the integrations over dz, dz' 
by the ones over dx, dx' ■ - The action for the propagation from t to t + e becomes 

ds^-i -V{q)^= l -^-eV^ + ix) , (10) 
where £ G [0, 1] ; S[q'} is calculated analogously. - Similarly, we evaluate the influence functional, 

rt+e rs 

!F[q, q', t + e] = T\q, q',t]+ ds du . . . 

Jt Jo 
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= J r [q,g',t]+ [ t+6 ds t du ... + 0(e 2 ) 
Jt Jo 

= ^[q,q',t}+0(e 2 ,e X ,ex') 

+e($ + X - $' - X') / du{v(t-u)[q(u)-q'(u)}+ir)(t-u)[q(u) + q'(u)}} , (11) 
Jo 

where we omitted the last two terms from Eq. (^). 

iv. Expanding consistently to 0(e 2 , ex, ex'); we obtain 

p($, t + e) ~ / dxdx'd^.d^ -L exp{^(x 2 - x' 2 )} (1 - zeV ($) + ieV($')) 

J Z7TC 2e 

P+**+** V((D(i exp {i{s[ q )-S[q'))}eM-^hQ',t)}p(^i,^0) 

J'\'..'\'.0 

■(l-e{$ + X-&-)d}£ du{u(t - u) [q{u) - q'{u)} +irj{t- u) [q{u) + q'{u)) }) ,(12) 

where a possibly unexpected normalization factor (27re) _1 appears on the r.h.s.; analogously to 
the factor (27re) -1 / 2 in the derivation of the Schrodinger equation from the path integral, where 
one employs only one intermediate coordinate integration, it is necessary here, in order that 
Eq. fljjD becomes an identity in the limit e — > . 

v. We observe that without the last factor under the remaining path integrals they would 
yield simply p($ + x, + x'?^) ■ Furthermore, the Gaussian integrals over d%, d%' effectively 
require x 2 >x' 2 — e - Therefore, properly expanding and collecting terms to 0(e), we obtain 
after performing these integrals the fully general result corresponding to the action of Eq. (Q): 

a tP (<s>, t) = -^[di - d£,]p{$, t) + - v(#)M$, *) 

_($_$') /" d $,d$: [*'*'* VqVq' exp{i{S[q] - S[j])}exp{-F[q,q',t]}p{$ i M,0) 

(-3 j 5 

• / du {v(t - M)k(w) - q'(u)] + in(t - u)[q(u) + q'(u)]} , (13) 
Jo 

where, of course, we cancelled the leading terms between the l.h.s. and the r.h.s. and divided 
the remaining ones by e, before taking the limit e — > . We recall that the quadratic potential 
terms have to be renormalized, 

i.e., one has to perform a mass shift here. 
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This completes our derivation of the master equation, which conforms with the general struc- 
ture of Eq. ([]]). Generally, we are still left with two formidable path integrals to be evaluated. 
However, depending on the actual form of the noise and dissipation kernels, the above result, 



Eq. dl3|) is a useful starting point for various approximation techniques, such as perturbation 



theory, mean field theory, variational techniques, or exact evaluations. 



The latter can be carried out in particular, if the action S is at most quadratic in $ |[T 
or else, if the kernels have zero range in time corresponding to the absence of non-Markovian 
memory effects. In the following section we will show how the previously studied Ohmic envi- 
ronment can be treated with the help of the above results. 

3 Example: Ohmic Environment 

In order to illustrate the explicit form of the master equation ([13]), we choose the simplest 
example of an Ohmic environment We begin with a quasi-continuous spectral density of 
the form UTTI 



I(u) = ^-u(-) k 1 exp{-u 2 /A 2 } , (15) 

7T \UjJ 

where A is a high-frequency cut-off, u another frequency scale, and 7 denotes the effective 
coupling constant between the $- field and the environment oscillators, cf. Eqs. (Q), (|5|). 

Then, considering the high-temperature limit T ^> A > uj and especially the Ohmic envi- 
ronment with k — 1 , we obtain the thermal noise and dissipation kernels according to Eqs. (£|) 

u(s) = 2 1 T5{s) , r)(s) = 7< J / (a) , (16) 

where the prime denotes the derivative of the 5-function in the last equation. Using these 
kernels in Eq. (|13"D, the path integrals simply reproduce p($, $',£) and the resulting equation 
becomes essentially local in time, 

{-\[di - dl\ + v(*) - v(&) - * 7 {2T[$ - $f + [$ - #][a, - a*]}} p($, t) 

+ 7 5(t) [$ - $'] J d^d^ [$< + Q'M^, 0) . (17) 
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We remark that herein a divergent term — 7^(0) [$ 2 — $' 2 ]p($, t) precisely cancelled the mass 
shift of Eq. (|14|); this has been the reason to introduce the counterterm J2n=i m n^n^ 2 /2 into 
the action, Eq. (|3|), in the first place. Furthermore, the singular initial value term in the last line 
of Eq. (17) arises similarly from evaluating by partial integration the last term proportional to 
7)(t — u) = j8'(t — u) in Eq. flllf ). This contribution may vanish for particular initial conditions. 
However, it has not been noticed in earlier work on the Ohmic environment by Caldeira and 
Leggett and by Unruh and Zurek as summarized in Ref. O]. It produces an initial 'jolt' 
which has been discussed in these studies of quantum decoherence. 

In order to cast the master equation (|TT| ) into the form of a sometimes more useful Fokker- 
Planck type transport equation, we introduce the Wigner transform 

W($, U,t) = J d0 exp{TI0}p($ + f , $ - f , t) . (18) 

With this definition and assuming that the potential V($) is a fourth order polynomial, we 
obtain the transport equation 

dtW($,n,t) = |n^-F , ($)9n + ^ ,,/ ($)9 n 3 + 47{9nn + T9 2 }|iy(<i>,n,t) 

-4-Kj5(t)5'(U) J d$' $'W« 0, 0) , (19) 

where the derivatives act on everything to the right and the primes on V(<&) denote differenti- 
ation w.r.t. $ . 

The various terms arising in Eq. ( |I9|) have the following significance: a) the l.h.s. together 
with the first two terms on the r.h.s. describe the classical Liouville flow of the phase space 
distribution W ($, II) ; b) the third term presents the quantum correction (for a non-polynomial 
potential one obtains here a well known sine series in powers of <9$<9n ); c) the fourth and fifth 
term, respectively, incorporate the previously obtained dissipative and diffusive contributions 
due to the (integrated out) Ohmic environment; d) the highly singular last term indeed induces 
an initial 'jolt' in the Wigner function, 

n, t = o+) = n, t = o_) - An^d'iu) J d$' $ / w / ($ / , o, t = o_) . (20) 

With this observation we conclude our (re-) derivation of more or less known features for the 
Ohmic environment and turn to a brief discussion of future applications of the method and 
main result presented in Sec. 2. 
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4 Discussion 

We derived the master equation for a homogeneous scalar field interacting with an environment 
composed of harmonic oscillators, which may alternatively be viewed as a free radiation field 
directly from the path integral representation of the density matrix propagator. In particular, 
we allowed for arbitrary self-interactions of the scalar field. We remark that our results apply 
equally well to a single nonrelativistic particle in an arbitrary external potential interacting 
with the environment, i.e. the case usually studied in quantum Brownian motion then, the 
variable $ is to be interpreted as the particle coordinate. 

As an immediate application of the results presently obtained we plan to study the phe- 
nomenon of semiquantum chaos in more realistic situations. 

Semiquantum chaos as described in Ref. |12[ (see also further references therein) arises 
in classically regular systems in the semiclassical regime. In particular, the time dependent 
Hartree-Fock approximation (or large-N expansion etc.) to the full quantum evolution, which 
in turn is regular, generally leads to nonlinearly coupled equations of motion for the one- and 
two-point functions. They have been shown to produce chaotic behavior under various circum- 
stances. 

Our aim is to show that the interaction with a suitably prepared environment, e.g. involving 
solid state quantum dots or QED microwave cavities, can drive the system, e.g. conduction 
band electrons or Rydberg atoms, into the semiclassical limit where semiquantum chaos is 
expected to manifest itself. 

As is obvious from the general result, Eq. (|TB|), or the more detailed calculation for the 
Ohmic environment which led to Eq. (|17|) and the Fokker-Planck type transport equation (19). 



these master equations all are linear in the density matrix or Wigner function. An important 
generalization of the approach of Sec. 2 would be to study the coupling between propagating 
fields. In this case we also expect nonlinear contributions to the effective 'one-body' master 
equation which are due to scattering terms. We hope to come back to the problem of propa- 
gating fields in the near future with a study of the entropy production due to bremsstrahlung 
emission during the initial stopping phase of high-energy heavy- ion collisions H. 
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